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ABSTRACT 

We use the observed polarisation properties of a sample of 26 powerful radio galaxies and 
radio-loud quasars to constrain the conditions in the Faraday screens local to the sources. 
We adopt the cosmological redshift, low-frequency radio luminosity and physical size of the 
large-scale radio structures as our ‘fundamental’ parameters. We find no correlation of the 
radio spectral index with any of the fundamental parameters. The observed rotation measure 
is also independent of these parameters suggesting that most of the Faraday rotation occurs 
in the Galactic foreground. The difference between the rotation measures of the two lobes 
of an individual source as well as the dispersion of the rotation measure show significant 
correlations with the source redshift, but not with the radio luminosity or source size. This is 
evidence that the small-scale structure observed in the rotation measure is caused by a Faraday 
screen local to the sources. The observed asymmetries between the lobes of our sources show 
no significant trends with each other or other source properties. Finally, we show that the 
commonly used model for the depolarisation of synchrotron radio emission by foreground 
Faraday screens is inconsistent with our observations. We apply alternative models to our 
data and show that they require a strong increase of the dispersion of the rotation measure 
inside the Faraday screens with cosmological redshift. Correcting our observations with these 
models for redshift effects, we find a strong correlation of the depolarisation measure with 
redshift and a significantly weaker correlation with radio luminosity. We did not find any 
(anti-)correlation of depolarisation measure with source size. All our results are consistent 
with a decrease in the order of the magnetic field structure of the Faraday screen local to the 
sources for increasing cosmological redshift. 
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1 INTRODUCTION 

The large-scale radio structure of radio galaxies and radio-loud 
quasars is created by the interaction of powerful jets with the 
gaseous environment of their host galaxies. The jets themselves 
are produced in the immediate vicinity of supermassive black holes 
at the centres of the host galaxies. Because radio radiation is not 
easily obscured, luminous radio-loud AGN are in principle ideal 
objects to study the co-evolution of supermassive black holes and 
their large-scale environments throughout the universe. However, 
determining the properties of the gas in the surroundings of these 
objects is not straightforward and a number of techniques have been 
applied in the past. 

Galaxy counts around individual host galaxies can be used to 
determine the galacti c environment, field, group or cl uster, of the 
radio source (e.g. Iffill & LillU 1 199 ll: I Wold et"ail l20001. Obviously 
the density of galaxies does not directly yield a constraint on the 
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properties of the gas in between the galaxies. A relation of galaxy 
to gas density must be assumed. 

The properties of the large-scale radio structure are shaped 
by the density of the gas the jets are ploughing into. The redshift 
distribution of, for e xample, the linear sizes of the lobes of radio 
galaxies of type FRII iFanaroff & Rilevtfl974t) in complete samples 
can be used to constrain the properties of the source environments 
fe.g. lKaiser & AlexandeJ.fl99^ :lBlundell et ailll999l) . Detailed fit¬ 
ting of two-dimensional radio maps can also be applied to individ¬ 
ual sources dkaise JS). However, the determination of the gas 
density from the source size is always model-dependent. Although 
most models for the evolution of radio lobes employ some form 
or other of ram pressure equilibrium between the jets and the re¬ 
ceding gas in the source environment, the model details introduce 
considerable uncertainty in the estimates of the gas density. 

The gas surrounding the lobes of radio galaxies consists of 
a magnetised plasma which acts as a Faraday screen for the ra¬ 
dio synchrotron emission of the lobes. By studying the polarisation 
properties of the radio emission, we can investigate the properties 
of this screen. The observable rotation measure caused by the Fara- 










2 Goodlet & Kaiser 


day screen depends only on the density of the plasma and also on 
the strength and orientation of the magnetic field in this material. 
Despite this degeneracy, we can hope to probe the properties of the 
gas in the source environments directly. The aim of the present pa¬ 
per is to use this tool to study a sample of 26 powerful radio sources 
defined in lGoodlet et alJ12004 hereafter G04). In order to break the 
degeneracy due to Malmquist bias between radio-luminosity and 
redshi ft observed in previous studies (e.g. iRronberg et all Il972t 
iLeahv et all f 1 9861) . we defined 3 subsamples of sources chosen 
from the 3C and 6C/7C catalogues. This selection enables us to 
investigate any trends of the polarisation properties of our sources 
with redshift and/or radio luminosity independently. 

In Section[2]we briefly describe our sample and the statistical 
methods we use. We also present the observational data that forms 
the basis for the subsequent analysis. Section[3]investigates trends 
of the polarisation properties with ‘fundamental' source parameters 
such as redshift, radio luminosity and linear size. Asymmetries be¬ 
tween the properties of the two lobes of each source are discussed 
in Section[4] In Section[5]we show that the polarisation properties 
of o ur sources are inconsistent with the commonly used model by 
ijBuni (Il966h for external Faraday screens. We then apply the more 
sophisticated models proposed bv iTribbie < 1199 1) to our data. Fi¬ 
nally, Section[6]summarises our main conclusions. 

Throughout this paper we use a cosmological model with 
Hq = 75 kms -1 Mpc -1 , £2 m = 0.35 and £2^ = 0.65. 


2 DATA AND STATISTICAL METHODS 
2.1 Sample description and observables 

Our sample of 26 sources is divided into three subsamples. Sample 
A wa s defined as a subsample chose n from the 6 CE lEal es et all 
1 19971) subregion o f the 6C survey lHales et all Il99cl) . and the 
7C III subsample Lacv et ah, 1999 1, drawn from the 7C survey 
IPoolev et al II 1 99? ). The selected sources have redshifts 0.8 < z < 
1.3, and radio luminosities at 151 MHz of 6.5 x 10 26 WHz~* < 
^151 MHz < 1-35 x 10 27 WHz _1 . Sample B was defined as a sub¬ 
sample from the revised 3CRR survey bv lLaing et alJ 0 983) con¬ 
taining sources within the same redshift range but with luminosities 
in the range 6.5 x 10 27 WHz-' < Pisimhz < 1-35 x 10 28 WHz-'. 
Sample C is also from the 3CRR catalogue; it has the same radio 
luminosity distribution as sample A, but with 0.3 < z < 0.5. For full 
details on the sample selection and the data reduction see G04. 

In the following we will use source properties derived from 
the observations: Spectral index, rotation measure and depolarisa¬ 
tion. All these properties were averaged over individual lobes. We 
define the spectral index, a, by S v v a . The depolarisation mea¬ 
sure, DM, is defined as the ratio of the percentage polarised flux at 
4.8 GHz and the percentage polarised flux at 1.4 GHz. The rotation 
measure, RM, is related to the degree of rotation of the polarisation 
position angle over a set frequency range, in our case from 1.4 GHz 
to 4.8 GHz. All source properties used here were determined for 
individual sources in G04. 

The angular resolution of our observations is limited by the 
need for information on the position angle of the plane of polari¬ 
sation at 1.4 GHz. For sources with lobes where the degree of po¬ 
larisation is uniform and the polarisation angle does not change on 
small scales, we can be confident that our measured depolarisation 
measures are accurate. However, some sources may exhibit polari¬ 
sation structure on small angular scales below our resolution limit 
and for these sources our measurements will be affected by beam 


depolarisation. Comparing our maps for the sources at low redshift 
presen ted in G04 with the higher resolution maps of [Gilbert et alJ 
t i2004i) . obtained mainly at 8 GHz, we find a few sources where 
beam depolarisation may play a role. Examples include 3C 351 and 
the northern lobe of 3C 46. Our general impression for this subsam¬ 
ple is that sources with distorted lobe structures may suffer from 
beam depolarisation. This should be borne in mind in the following 
analysis. Unfortunately, in the absence of better data, particularly 
for the sources at higher redshifts, it is not possible to quantify this 
effect. 

The difference in the spectral index, depolarisation and rota¬ 
tion measure over the lobes is given by da, d DM and ARM respec¬ 
tively. All differential properties are taken to be the difference of 
the averages over individual lobes. The rms-variation or dispersion 
of the rotation measure across a lobe is defined by Orm- 

In G04 we argue that the average rotation measure of the 
sources, RM, is mainly caused by the interstellar medium in the 
Galaxy rather than the material in the vicinity of our sources. How¬ 
ever, fluctuations in the Faraday screen in the source environment 
will give rise to the observed variations of RM across the source on 
smaller spatial scales as measured by ARM and Crm, as well as the 
depolarisation, DM. Thus we need to correct for the redshift depen¬ 
dence of these quantities. We correct ARM and Crm to the sources’ 
frame of reference by multiplying the measured values by (1 + z) 2 . 
Note that this factor only corrects for the wavelength dependence 
of RM. It does not correct for the resolution effects discussed in 
Section|5] We also adjust DM and ADM to their values appropriate 
for a uniform cosmological redshift z = 1. The latter adjustment re¬ 
quires a model for the structure of the Faraday screen local to the 
source which gives rise to these meas urements. We initially apply 
the commonly used model presented in lBur 3 il966t hereafter B66), 
but discuss alternative, potentially more accurate models that also 
affect our measurements of CJrmz in Section[5] The adjusted quan¬ 
tities are labelled ARM Z , ctrmz, DM : and ADM Z . We summarise the 
observed properties of all sources in our sample in TableQ 
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Name 

z 

P 151 / WHz-‘ 

ID 

D / kpc 

S 4.8 / mJy 

a RM / rad m 2 

Ormz / rad m - 

DM Z 

6C 0943+39 

1.04 

1.0 

X 

10 27 

w 

56 

31.0 

-0.72 

1.6 

67.0 

2.36 






E 

45 

42.0 

-1.01 

-19.1 

422.0 

4.38 

6C 1011+36 

1.04 

1.1 

X 

10 27 

N 

175 

44.3 

-0.88 

30.8 

53.3 

1.60 






s 

220 

18.7 

-0.89 

12.6 

42.0 

1.07 

6C 1018+37 

0.81 

8.8 

X 

10 26 

NE 

279 

46.4 

-0.85 

0.85 

20.0 

1.15 






sw 

294 

28.7 

-0.84 

11.2 

9.5 

1.82 

6C 1129+37 

1.06 

1.1 

X 

10 27 

NW 

63 

46.5 

-0.85 

-19.3 

229.2 

2.86 






SE 

69 

73.2 

-0.90 

0.0 

93.4 

6.35 

6C 1256+36 

1.07 

1.3 

X 

10 27 

NE 

69 

57.8 

-0.79 

5.9 

93.4 

1.42 






SW 

62 

101.4 

-0.87 

15.4 

60.0 

1.01 

6C 1257+36 

1.00 

1.1 

X 

10 27 

NW 

134 

43.5 

-0.71 

-115.3 

40.0 

1.27 






SE 

150 

20.5 

-1.06 

-115.6 

64.0 

1.90 

7C 1745+642 

1.23 

9.5 

X 

10 26 

N 

66 

23.5 

-0.86 

— 

— 

2.71 






S 

99 

33.8 

-0.91 

12.8 

15.4 

1.00 

7C 1801+690 

1.27 

1.0 

X 

10 27 

N 

89 

8.8 

-0.97 

44.8 

82.5 

1.19 






S 

113 

28.7 

-0.81 

20.8 

51.5 

1.73 

7C 1813+684 

1.03 

7.1 

X 

10 26 

NE 

282 

15.0 

-0.92 

13.9 

350.3 

1.02 






SW 

152 

30.2 

-0.79 

-68.4 

171.8 

1.15 

3C65 

1.18 

1.0 

X 

10 28 

W 

83 

524.0 

-1.00 

-82.6 

101.7 

6.12 






E 

68 

240.9 

-1.03 

-86.1 

46.1 

1.42 

3C68.1 

1.24 

1.1 

X 

10 28 

N 

178 

667.2 

-0.82 

-26.6 

139.0 

1.45 






S 

223 

36.8 

-1.08 

57.9 

300.1 

2.69 

3C252 

1.11 

7.2 

X 

10 27 

NW 

159 

178.7 

-1.00 

15.7 

89.9 

1.10 






SE 

257 

80.0 

-1.10 

58.5 

202.1 

2.50 

3C265 

0.81 

6.5 

X 

10 27 

NW 

331 

224.0 

-0.62 

42.2 

71.8 

1.52 






SE 

216 

318.9 

-0.78 

32.8 

74.7 

1.35 

3C 268.1 

0.97 

1.0 

X 

10 28 

E 

151 

262.3 

-0.92 

21.7 

41.9 

1.52 






W 

194 

2296.6 

-0.58 

26.8 

51.6 

1.82 

3C267 

1.14 

1.0 

X 

10 28 

E 

158 

184.0 

-1.17 

-9.6 

111.7 

1.43 






W 

129 

479.2 

-0.83 

-21.5 

107.2 

1.08 

3C280 

1.00 

1.2 

X 

10 28 

E 

51 

326.0 

-1.01 

-37.7 

80.4 

1.83 






W 

61 

1289.2 

-0.76 

-7.5 

6.0 

1.54 

3C324 

1.21 

1.3 

X 

10 28 

NE 

50 

432.6 

-1.05 

22.1 

112.3 

2.15 






SW 

52 

166.0 

-1.14 

43.0 

82.0 

2.75 

4C 16.49 

1.29 

9.9 

X 

10 27 

N 

89 

120.0 

-0.88 

-4.3 

295.8 

1.10 






S 

51 

142.3 

-1.32 

0.9 

299.8 

1.16 

3C16 

0.41 

7.4 

X 

10 26 

SW 

200 

484.9 

-0.97 

-4.3 

34.6 

1.16 






NE 

187 

22.1 

-0.87 

— 

— 

1.22 

3C42 

0.40 

7.5 

X 

10 26 

NW 

87 

353.9 

-0.87 

-2.4 

27.6 

1.01 






SE 

92 

450.6 

-0.86 

5.0 

27.8 

1.00 

3C46 

0.44 

7.9 

X 

10 27 

NE 

512 

162.7 

-0.94 

-4.8 

8.1 

1.32 






SW 

376 

173.7 

-0.91 

-2.9 

8.7 

1.17 

3C299 

0.37 

6.9 

X 

10 26 

NE 

71 

876.5 

-0.93 

-126.3 

142.5 

1.15 






SW 

36 

53.7 

-0.71 

16.0 

12.9 

1.04 

3C341 

0.45 

8.8 

X 

10 26 

NE 

280 

123.8 

-0.85 

20.3 

21.0 

1.11 






SW 

145 

265.9 

-1.00 

18.2 

25.2 

1.13 

3C351 

0.37 

7.6 

X 

10 26 

NE 

158 

1093.9 

-0.80 

1.0 

23.5 

1.02 






SW 

188 

77.1 

-0.93 

4.4 

20.3 

1.24 

3C457 

0.43 

9.6 

X 

10 27 

NE 

568 

208.6 

-1.02 

— 

— 

1.00 






SW 

519 

290.5 

-0.94 

— 

— 

1.05 

4C 14.27 

0.39 

8.8 

X 

10 26 

NW 

103 

107.0 

-1.06 

-13.0 

10.4 

1.09 






SE 

76 

124.8 

-1.17 

-17.3 

9.3 

1.01 


Table 1 . Observational properties of the sources in our sample. P 151 is the radio luminosity at 151 MHz, D the physical size of individual lobes and S 4.8 is the 
radio flux at 4.8 GHz. Ormz is the dispersion of the rotation measure corrected for redshift effects. DM Z is the depolarisation adjusted to a common redshift of 
Z — 1 using the model of|Burrilil966). All data taken from G04. 
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2.2 Statistical methods 

In this paper we use the Spearman Rank test (SR), the Partial Spear¬ 
man Rank test (PSR) and Principal Component Analysis (PCA). 
SR is a non-parametric correlation test which assigns a rank to 
given source properties, say X and Y, and computes the coef¬ 
ficient of the correlation, rxy, between the ranks. By definition 
— 1 < rxY < 1, where negative values of rx y indicate an anti¬ 
correlation. The null hypothesis that no correlation is present corre¬ 
sponds to rxY = 0 an d in general rx y has a Student-t distribution. 
A (anti-)correlation with a confidence level exceeding 95% will re¬ 
sult in a Studen t-t valu e greater than 2 for the size of our sample 
(e.g. lRilev et all 1997"). 

PSR tests whether a correlation between properties X and Y is 
due to independentmdividual correlations of X and Y with a third 
property Z (e.g. lMacklinlll983) . The correlation coefficient, rxy,z> 
increases (decreases) from zero if the correlation (anti-correlation) 
is predominantly between X and Y, independent of Z. The positive 
significance of the (anti-)correlation is measured by Dxy z- 

PCA is a multi-va riate statistical test (e.g. iDeemin e,® 
lEfstathiou & Falill984l) . PCA is a linear self-orthogonal transfor¬ 
mation from an original set of n objects and m attributes forming 
a nx m matrix, with zero mean and unit variance, to a new set of 
parameters, known as principal components. The principal com¬ 
ponents of the new dataset are all independent of each other and 
hence orthogonal. The first component describes the largest varia¬ 
tion in the data and has, by definition, the largest eigenvalue of the 
n x m matrix. In physical terms the magnitude of the eigenvalue 
determines what fraction of the variance in the data any correla¬ 
tion describes. Thus a strong correlation will be present in the first 
eigenvector and will only strongly reverse in the last. PCA in its 
simplest terms is an eigenvector-eigenvalue problem on a trans¬ 
formed, diagonalised and standardised set of variables, n objects 
will create an eigenvector-eigenvalue problem in n — 1 dimensions. 
To avoid problems with the interpretation of our results, we do not 
use more than 4 parameters in any part of the analysis. 


3 TRENDS WITH REDSHIFT, LUMINOSITY AND SIZE 

The aim of our analysis is to investigate any trends with redshift 
and/or radio luminosity of observables that may depend on the 
properties of the gas surrounding the radio lobes of our sample 
sources. Another important parameter that may influence our re¬ 
sults is the physical size of the radio lobes. Because we expect 
significant gradients for the gas density in the typical environ¬ 
ments of sources, the properties of small sources will clearly be 
shaped by very different conditions than those of large sources, 
even if the overall gaseous environment is the same. Many stud¬ 
ies have investigated the dependence of the average size of sources 


pies of sources (e.g. Singal, 1988 

; Barthel & Miley, 1988; Kanahi 

1989; Singail993 iNeeser et al„ 

1995; Kaiser & Alexander, 1999 

Blundell et al. 1999). This is not the purpose of the study presented 


here. Our sample is not complete and sources were selected in order 
to cover an equal range of physical sizes within our sample limits 
of redshifts and radio luminosities. 

In the following we will refer to the redshift, z, radio luminos¬ 
ity at 151 MHz, P, and the physical size, D, as the ‘fundamental’ 
source parameters. In order to better understand the correlations of 
other observables with the fundamental parameters, we first need to 
investigate any correlations between them. This analysis also serves 
as a useful test of the robustness of our sample selection. 


X,Y 

rxY 

Student-t 

Z 

rx Y,z 

Dx Y,Z 

z,P 

0.70 

4.8 

D 

0.69 

4.0 

Z,D 

-0.23 

-1.2 

P 

0.04 

0.20 

P.D 

-0.13 

-0.64 

z 

-0.20 

0.93 


Table 2. Spearman Rank coefficients for correlations between the funda¬ 
mental parameters. 


X, Y 

rx Y 

Student-t 

X, Y 

rx Y 

Student-t 

ot, z 

-0.13 

-0.62 

da, z 

0.15 

0.77 

a, P 

0.04 

0.19 

da, P 

0.27 

1.35 

a. D 

0.21 

1.05 

da, D 

-0.21 

-1.05 


Table 3. Spearman Rank coefficients for the correlations between spectral 
index, a, and difference of spectral index between lobes, da, and the fun¬ 
damental parameters. 


Table |2| illustrates the relations between the fundamental pa¬ 
rameters in our sample. Clearly, redshift and radio luminosity are 
strongly correlated, but this correlation simply follows from our 
sample selection. Using the statistical techniques described above, 
we can isolate this correlation from others. In fact, it provides us 
with a ‘bench mark’ for other correlations. 

The SR results also imply a weak anti-correlation between the 
redshift of a source and its physical size. However, the PSR results 
show that the relation between the three fundamental parameters 
is more complex than a simple anti-correlation between z and D. 
We also found that the anti-correlation disappears if the two ex¬ 
ceptionally large sources 3C 46 and 3C 457, both at low redshifts, 
are removed from the sample. In practice, the z-D anti-correlation 
is too weak to introduce any significant bias into our results. We 
therefore conclude that our sample contains a sufficient range of 
physical sizes at each redshift and radio luminosity. 

In the following we investigate not only the observed source 
properties, but also differences of these properties between the two 
lobes of each source. Clearly there are twice as many direct mea¬ 
surements of, say, the spectral index of individual lobes, than mea¬ 
surements of the difference of the spectral index between the lobes. 
To avoid complications with the different number of measurements, 
we always use values for all source properties averaged over the en¬ 
tire source. 


3.1 Spectral index 

The dependence of the spectral index, a, on the fundamental pa¬ 
rameters for radio galaxie s of type F RII has been inves tigated 
by various authors. IVeron et alJ C 1972i ) and (o nuor 1 il989i) find a 
strong spectral index-radio luminosity correlation, but no corre- 
spondi ng correlation with redshift. Interestinelv lAthreva & Kanahil 
1 1999 ) find a strong correlation of spectral index with red¬ 
shift, but no si gnifica nt correlation with radio luminosity. Finally, 
{Blundell et aL 1199 9) find that a is anti-correlated not only with 
radio-luminosity, but also with source size. They found that in gen¬ 
eral, larger sources had a steeper spectral index. 

In our sample there is no clear trend of the spectral index with 
any of the fundamental parameters (Tables[3]and[4j. The difference 
in the spectral index of the lobes, da, also shows no significant 
correlations with the fundamental parameters. The absence of the 
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Eigenvector 



i 

2 

3 

4 

z 

0.64 

-0.25 

0.06 

-0.73 

p 

0.58 

-0.16 

0.53 

0.60 

D 

-0.47 

-0.12 

0.82 

-0.31 

a 

-0.21 

-0.95 

-0.21 

0.13 

Variation 

45% 

25% 

21 % 

10 % 

z 

0.56 

-0.08 

0.54 

-0.62 

p 

0.55 

0.44 

0.26 

0.66 

D 

-0.40 

0.84 

0.19 

-0.30 

da 

0.48 

0.30 

-0.77 

-0.29 

Variation 

51% 

21 % 

18% 

10 % 


Table 4. PCA results for the spectral index and the fundamental parameters. 


X, Y 

OCY 

Student-t 

X, Y 

rxY 

Student-t 

RM, z 

0.17 

0.79 

ARM, z 

0.50 

2.65 

RM, P 

0.27 

1.31 

ARM. P 

0.23 

1.06 

RM, D 

-0.05 

-0.23 

ARM, D 

-0.22 

-1.05 

CRMz. Z 

0.68 

4.24 




ORMz, P 

0.35 

1.71 




ORMz, D 

-0.30 

-1.46 





Table 5. Spearman Rank coefficients for the correlations between rotation 
measure, RM, the difference of the rotation measure between lobes, d RM, 
and the dispersion of the rotation measure, <3rm z . and the fundamental pa¬ 
rameters. 

trends claimed in the literature from our sample is probably due to 
the small sample size. 

3.2 Rotation measure 

The SR results in Table[5]show that there are no strong correlations 
of the fundamental parameters with rotation measure, RM. In every 
case the correlation is below 95% significant. The PCA results in 
Table [6] also demonstrate that the RM of a source does not depend 
on the size, redshift or the radio luminosity of the source. This is 
consistent with our interpretation that the observed RM is caused 
by the Galactic magnetic field and not by Faraday screens local to 
the sources lC104: lLeah\l 1 9871) . 

The difference in the rotation measure, d RM : , between the two 
lobes of a given source and the dispersion of the rotation measure, 
Ormz) are probably caused locally to the source and have therefore 
been corrected to the sources’ frame of reference (G04). There is 
a ARM z -z correlation present in the SR results, (see Table [5j, but 
it is weaker in the PCA results where the correlation is strongly 
reversed in the 3rd eigenvector (see Table|6j. FigureQshows that in 
principal we should find a strong d RM z -z correlation. However, the 
presence of the source 3C 299 at low redshift, but with a large value 
of d RM Z weakens this trend. The large value of d RM Z for 3C 299 
is due to a lack of polarisation information in one lobe. Removing 
this source strengthens the correlation giving r ( \rm, z = 0.632 which 
corresponds to a significance of 99.98%. 

°RMz measures variations of the Faraday screen on scales 
smaller than d RM Z , which may in principle be influenced by the 
Galactic Faraday screen. However, it is interesting to note that 
d RM Z correlates strongly with ctrmz ( see Figure^. So the two pa- 
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0.03 

ORMz 

0.52 

-0.06 

0.71 

-0.47 

Variation 

51% 

22 % 

19% 

8 % 


Table 6. PCA results for the rotation measure and the fundamental param¬ 
eters. 
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Figure 1. Redshift against difference in rotation measure, d RM Z across 


lobes of a source, d RM- has been shifted to the sources’ frame of reference. 


rameters may sample the same, local Faraday screen even if d RM Z 
may still contain some Galactic contribution. 

(Trmz shows a strong correlation with redshift, > 99.9% sig¬ 
nificance, and a weak anti-correlation with D. The Grmz-z cor¬ 
relation also strengthens with the removal of 3C299, becoming 
r a rmz,z = 0-767 which corresponds to a significance exceeding 
99.99%. These correlations are more obvious when the PCA results 
are considered (see Table|6|. The first eigenvector contains 51% of 
the variation in the data and shows a strong z-Ormz correlation and 
a weaker D-Ormz anti-correlation. The D-Ormz anti-correlation re¬ 
verses in the third eigenvector indicating that it is a considerably 
weaker trend compared to the z-Grm correlation which only re¬ 
verses in the last eigenvector. The PSR results in Table ^confirm 
that there are no significant correlations with radio luminosity for 
both d RM Z and (Trmz-_ 

In a survey of 27 high redshift sources lPentericci et alj 1 2001ill 
found that the Faraday rotation was independent of size or ra- 
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Table 8. Spearman Rank coefficients for the correlations between depolar¬ 
isation, DM Z and the difference of the depolarisation between lobes, d DM Z , 


and the fundamental parameters. 


Table 7. PSR coefficients for correlations of d RM Z and Ormz with the fun¬ 
damental parameters. 
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0.75 
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52% 

22 % 
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Table 9. PC A results for the depolarisation and the fundamental parameters. 


Figure 2. The dispersion of the rotation measure, (Jrm z , against the differ¬ 
ence in rotation measure across the lobes of a source, d RM Z . 

dio luminosity, but also showed that the number of sources with 
high level s of Faraday rot ation increased with redshift. The re¬ 
sults from IPentericci etlrt] i2000l) use the RM of a source in the 
sources’ frame of reference. As noted in G04 and above, the 
RM observed in our sources is dominated by contributions from 
the Galaxy. Theref ore we cannot d irectly compare our results 
wi th those of | Pen tericci et"aD fcOOOl) . Also note that the sample 
of IPentericci et alJ feOOCt) contains only sources at considerably 
higher redshifts than our s ample . 

The trend reported in lPentericci et a fl Eooa indicates that the 
strength of the magnetic field and/or the density of the environment 
is increasing with redshift. Which effect dominates cannot be de¬ 
cided from RM measurements alone. Similarly, we use Ormz and 
d RM Z which also do not provide us with a direct measurement of 
the density or the strength of the magnetic field in the vicinity of 
our sources. However, our indicators depend strongly on the de¬ 
gree of order in the direction of the magnetic field in the local Fara¬ 
day screen. Thus our results point towards a trend for more chaotic 
magnetic field structures in the source environments at higher red¬ 
shift. 

3.3 Depolarisation 

Depolarisation arises from structure of the magnetic field in the 
Faraday screen at scales below the telescope resolution. In the pre¬ 
vious Section we already found that the measured dispersion of 
RM is correlated with the redshift of the source, indicating that the 
small-scale structure in the RM arises from a Faraday screen lo¬ 
cal to our sources. Therefore we expect that the depolarisation also 
correlates with redshift in our sample. Clearly we must also correct 
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3.25 
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2.77 

P, d DM Z 
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0.05 

0.23 


Table 10. PSR coefficients for correlations of DM Z and d DM- with the fun¬ 
damental parameters. 


12 3 4 5 

DM, 

Figure 3. Redshift against the average depolarisation of a source. The de¬ 
polarisation has been shifted to a common redshift. z = 1 ■ 
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Figure 4. Radio-luminosity against the average depolarisation. The depo¬ 
larisation has been shifted to a common redshift, z = 1 ■ 


our measured DM values for redshift effects. In this Section we ap¬ 
ply the commonly used model of B66 to adjust DM and d DM to a 
common redshift of z = 1 for all our sources to allow for a direct 
comparison. In Section[5]we will investigate alternative models. 

There is little difference between the redshift and the radio 
luminosity correlations with the average adjusted depolarisation, 
taken over the source as a whole, DM Z , (see Table|8j as both cor¬ 
relations show significance levels in excess of 99.9%. There is also 
an indication of a DM Z -D anti-correlation, but this trend is not sig¬ 
nificant. We use the Partial Spearman Rank test to determine if the 
P-DM Z correlation can be explained by the independent z-DM z and 
z-P correlations (see Table llOt . Again we find little difference be¬ 
tween the redshift and radio-luminosity correlations. This result 
is also confirmed by the PC A results in Table [9] Figure [3] shows 
that the higher redshift sources display, on average, a larger degree 
of depolarisation and also a much larger spread of depolarisation 
at any given redshift compared to their low redshift counterparts. 
There is little difference in the range of DM Z between sources with 
very different radio luminosity (see Figure[4j. However, this result 
may be caused by only two sources of low radio luminosity and 
exceptionally high values of DM Z . Previous studies of the correla¬ 
tion of depolarisation with redshift and/or radio luminosity are also 
inconlusive. lMorris & Tabari j 1 9731 ) found depolarisatio n to corre¬ 
late with radio luminosity whereas TKronbere et alJll972lt found the 
correlation was predominantly with redshift. In view of the stronger 
correlation of Grmz with redshift it is surprising that DM Z does not 
show a correlation with redshift that is significantly stronger than 
the correlation with radio luminosity. However, we show in Sec- 
tion[5]that this result depends on which model for the local Faraday 
screen is used to correct for redshift effects. 

The difference in the depolarisation over the source. dDM z , 
shows a stronger trend with redshift than the overall depolarisa¬ 
tion, DM Z (Tables [8] and The PSR results, Table fTol indicate 
a significant correlation between redshift and d DM Z at any given 
radio luminosity. The corresponding P-dDM z correlation is consid¬ 
erably weaker at a given redshift. 

Table [8] shows a weak DM Z -D anti-correlation. The PCA 
results in Table [9] confirm that this result is not very signifi- 
cant. Howev e r, a DM 7 -D anti - correl a tion has been reported by 
IStroml ll973^JStron^JagM^ (j_98jj); |Pedelty_etaI] (j_989 l here - 
after P89l. I B es‘HiHin"^^ 99n and Ilshwara- Chandra et alJ 1 1 99~ij) . 
but was not found bv IPennett-ThorneT^^fifl . The observations by 


[Strom (1973 ). IStrom & Jaeersl 1 1988h and Ifshwara-Chandra et alJ 
( 1998 ) are at a lower resolution than our measurements and thus 
their different findin gs may be_due to this resolution difference. The 
observations bv lBest et al. ill 999 ) are at a higher resolution, but also 
at higher frequencies (4.8 GHz to 8.4 GHz) and so we would expect 
any depolarisatio n trends to be stronger in this sample. Both P89 
and iDennett-Thornd (1996) however observe the same frequency 
range and with the same spatial resolution as in our sample. Even 
when our results from sample B, which are taken at the most simi¬ 
lar redshift and radio-luminosity range to those of P89. are analysed 
separately, we still find no trend with size. There is no obvious dif¬ 
ference in the selection of the P89 sample and our sample B, except 
that P89 chose their sample from sources which had strong emis¬ 
sion lines. By preferentially selecting sources with strong emission 
lines P89 have chosen sources with jet axes close to our line of 
sight. This selection may strengthen any D-DM correlation found. 
However, we do not have a convincing physical explanation why 
our results differ from the ones of P89 because of this selection 
criterion. 

It is usually assumed that radio sources are located in stratified 
atmospheres. Therefore a small source will be embedded in denser 
gas which acts as a more efficient Faraday screen. As the source 
expands the lobes will extend beyond the denser inner atmosphere, 
thus reducing the amount of depolarisation observed. In the begin¬ 
ning of Section 0 we showed that there is no significant z-D anti¬ 
correlation in our sample. Therefore the lack of a pronounced DM Z - 
D anti-correlation suggests that there is no significant difference 
in the source environments at any redshift as the sources become 
larger. This may be evidence that the environments of our sources 
are relatively homogeneous with little stratification. The lack of an 
anti-correlation between the source size, D, and d DM Z (see Tables 
[8]and|9} is consistent with this idea. 


4 SOURCE ASYMMETRIES 

In Section 0 we investigated the trends of observables averaged 
over the entire structure of our sources with respect to redshift, to¬ 
tal low-frequency radio luminosity and overall source size. In this 
Section we study in detail the asymmetries between the lobes of 
our sources, in particular how differences in the lengths and ra¬ 
dio fluxes of lobes of individual sources relate to asymmetries in 
other observable properties. In the following discussion on possi¬ 
ble mechanisms that could give rise to asymmetries between the 
two lobes of a source, we assume that the two jets of a given source 
are not significantly different from each other in terms of their en¬ 
ergy transport rate. 

The angle at which a source is orientated to the line-of-sight 
affects the projected lobe length. Because of finite light travel times 
sources orientated at a small angle will have lobes that appear more 
asymmetrical than sources at a large angle, assuming there are no 
differences between the environments of the two lobes. Assuming 
a source is oriented at an angle 0 to the lin e-of-si ght, then the ratio 
of lobe lengths, D 1 /D 2 , is given by (e.g. Longair & Ri ievHT97^) 

D\ _ 1 + (vo/c) cos0 
Di 1 — (v 0 /c) cos0’ 

where v 0 is the velocity at which the hotspot is moving away from 
the nucleus and c is the speed of light. The orientation of a source 
will also induce asymmet ries in the the rotation measure and the 
depolarisation (e.g. lGarrington & Conwavill99ll ). In a simple ori¬ 
entation model emission from the lobe pointing away from the 
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observer will have a longer path length through the local Fara¬ 
day screen and thus any depolarisation measurements of this lobe 
would be larger than the lobe pointing towards the observer, which 
is referred to as the Laing-Garrington effect. 

The motion of material in the radio hotspots at the end of the 
lobes of FRII-type sources is likely to be relativistic. Relativistic 
beaming can therefore also cause asymmetries in a source, inde¬ 
pendent of path length through the Faraday screen. As the hotspot 
becomes increasingly more beamed, i.e. 0 —> 0, the lobe ratio in¬ 
creases. The flat radio spectrum of the hotspot will begin to dom¬ 
inate the flux and spectral index in the beamed lobe. Thus as rela¬ 
tivistic beaming becomes more dominant in a source, the beamed 
lobe will become brighter, longer and its spectrum flatter than the 
receding lobe. iLiu & Poolevl i 199la ) also find for sources in their 
sample that the least depolarised lobe has a flatter spectrum. 

Finally, differences in the environments of the two lobes may 
also introduce asymmetries. For example, theoretical models for 
the evolution of radio lobes predict that denser environments lead 
to shorter and brighter lobes with steeper spectral indices (e.g. 
iKaiser et all ll997 ). A jet with a given energy transport rate will 
expand more slowly in a dense environment compared to the situ¬ 
ation in less dense surroundings. The pressure in the lobe as well 
as the strength of the magnetic field will be higher. The increased 
strength of the magnetic field will lead to increased energy losses 
of the radiating relativistic electrons, resulting in a steeper spectral 
index. Thus a source situated in an environment with an asymmet¬ 
rical density distribution will develop asymmetrical lobes. In the 
following we will use the lobe asymmetries observed in our sample 
in an attempt to distinguish between these possibilities. 

In Section we only used the overall physical size of our 
sources, D. We now need to measure the lengths of individual lobes. 
In about half of our sources we do not detect a radio core at our two 
observing frequencies. In these cases, wherever possible, we used 
previously published maps to establish the core position. However, 
in several sources, most notably 3C 16, no core is reported in the lit¬ 
erature. We estimated the location of the core position by using the 
most likely core position from the literature (Liu & Poolevt l 1991 b). 
This approach is not ideal, but as our maps are not of very high 
resolution, the values calculated for the lobe length will only be 
approximations to the true lengths in any case. 

In Section [3] we adjusted the measured depolarisation to a 
common redshift of z = 1. For this correction the we had to as¬ 
sume a specific underlying structure of the Faraday screen. Here 
we compare the two lobes of a given source and so we can use 
the uncorrected depolarisation to avoid any model dependence. 
This choice is consistent with previ ous studies fP89 lLiu & Pooled 
ll991al : lGarrington & Conwavll99lh which also used the observed, 
uncorrected depolarisation. 


4.1 Lobe length and radio flux 

Figure [5] shows the ratio of lobe lengths as a function of overall 
physical size of the sources. Most of our sources are quite symmet¬ 
ric and there is no tendency for smaller sources to show larger lobe 
length asymmetries than larger sources. 

We use the radio flux at 4.8 GHz in our study of the source 
asymmetries, since it has lower noise levels than the 1.4 GHz data 
(G04). Figure [6] demonstrates that for all sources the flux of the 
brighter lobe is tightly correlated with the flux of the fainter lobe. In 
a few cases the flux of the fainter lobe is much smaller than the flux 
of the brighter lobe. Three of these very asymmetrical sources are 
quasars, 3C 351, 3C 68.1 and 3C 16, the jet axes of which are pre- 
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Figure 5. Plot of the lobe ratio against the physical size of the entire radio 
structure of our sources. 




Figure 6. Plot of the flux of the fainter lobe against the brighter lobe at 
4.8 GHz. The line shows the location of sources for which the luminosities 
of the two lobes are identical. 
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Figure 7. Plot of the lobe ratio against the flux ratio at 4.8 GHz, omitting 
the sources with an exceptional flux asymmetry, 3C299, 3C 268.1, 3C 16, 
3C351 and 3C68.1. 
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Figure 8. Plot of the spectral index ratio against the lobe ratio. 


sumably closely aligned with our line of sight probably leading to 
the large flux asymmetry. The radio galaxies 3C299 and 3C 268.1 
also show a large flux as ymmetry. The pr onounced asymmetry of 
3C 299 was also noted bv lMcCarthv et alJll995l) . who suggest that 
despite its low redshift the source shows characteristics of a high 
redshift (z >1) radio galaxy. We also find 3C299 to have a Grm 
comparable to the high redshift samples. 

lHarvanek & Stock3 i20Q2h noted that 3C 299 has lobes of very 
different lengths. In fact, it is the source with the largest length 
asymmetry in our sample. The other four sources showing a large 
flux asymmetry do not have unusual lobe ratios. However, it should 
be noted that the cores in 3C 68.1 and 3C 16 have no t been detected 
in our observations or previously. iBridle et a] J j 1 9941) do marginally 
detect a core in deeper observations of 3C68.1 and we have used 
their core position in determining the lobe lengths of this source. 
To further aggravate the problem, both, 3C 16 and 3C68.1, have 
very faint lobes and so it is hard to determine where exactly the 
lobe terminates. However, it is worth noting that when the 5 un¬ 
usual sources are removed from consideration there is no signifi¬ 
cant trend for the sources with a large flux ratio to also show strong 
length asymmetries (see Figure^. This argues against the interpre¬ 
tation that the source asymmetries are caused by differences in the 
gaseous environments of the two lobes of a source. 

We do not find any trends of the flux asymmetries with either 
redshift or radio luminosity. 

4.2 Spectral index 

The spectral index ratio shows no correlation with the lobe length 
ratio (see Figure|8j. If differences in the gaseous environment of the 
lobes were the source of the observed asymmetries, then we would 
expect that the lobe located in the denser environment is shorter and 
has a steeper radio spectrum fe.g. lKaiser et all 11997 1. This result 
therefore argues against different lobe environments. However, it 
is consistent with relativistic beaming of the flat-spectrum hotspot 
emission, if the lobe length asymmetry is due to only small differ¬ 
ences in the gas density of their environments. 

Neglecting the five sources with unusually large flux asym¬ 
metries, we find that in sources with a large flux asymmetry the 
brighter lobe has a flatter spectral index than the fainter lobe (see 
Figure |9j. For more symmetric sources the range of the spectral 
index ratio is larger with a weak tendency for the brighter lobe to 
have a steeper spectral index than the fainter lobe. Both results are 


Figure 9. Plot of the spectral index ratio against the flux ratio at 4.8 GHz. 
omitting the sources with an exceptional flux asymmetry, 3C 299, 3C 268.1, 
3C16, 3C351 and 3C68.1. 
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Figure 10. Plot of the depolarisation measure ratio against the spectral 
index ratio . S ources above the dotted line follow the trend proposed by 
iLiu & Poolev l ( 1991a ). 


again consistent with asymmetries arising from relativistic beaming 
of the hotspot emission. Lobes for which beaming is more impor¬ 
tant will be dominated by their bright hotspots which have a flat 
radio spectrum. Relativistic beaming is less important for sources 
oriented close to the plane of the sky. For these objects small differ¬ 
ences in the density of the lobe environments may lead to a brighter 
lobe with a steeper spec trum located in the denser surroundings. 

ILiu & Poolevl 1 199la ) study a sample of 13 sources, for 12 of 
which the lobe with the steeper spectrum is also the most depo¬ 
larised lobe. Figure fTol illustrates that only 14 of our 26 sources are 
consistent with this trend. However, our findings are consistent with 
those of Ilshwara-Chandra et al .12001) who also find only 58% of 
radio galaxies and 59% of quasars in their much larger sample to 
agree with the result ot lLhi & Poolevi (1991a). The distribution in 
Figure flol suggests that sources with a larger spectral index asym¬ 
metry are more likely to follow the proposed trend. Again this is 
consistent with the asymmetries arising from relativistic beaming 
of the hotspot emission. In this scenario, the lobe with the steeper 
spectral index points preferentially away from the observer. The 
radio emission from this lobe has to traverse a longer path through 
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Figure 11. Plot of the depolarisation measure ratio against the lobe length 
ratio. For sources above the dotted line the larger lobe is also more depo¬ 
larised. 
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Figure 12. Plot of the depolarisation measure ratio against the flux ratio at 
4.8 GHz, omitting the sources with an exceptional flux asymmetry, 3C299, 
3C268.1, 3C 16, 3C351 and 3C68.1. For sources above the dotted line the 
brighter lobe is also more depolarised. 


the surrounding medium which acts as a Faraday screen, leading to 
greater depolarisation. 


4.3 Depolarisation 

In Figure fTTI we compare the depolarisation measure ratio with the 
ratio of the lobe lengths. For 16 of our 26 sources the larger lobe is 
also more depolarised. Our result disagrees with the findings of P89 
andj Liu & Poolevl j 1991 all , but is consistent with iDennett-Thornel 
(1996). The cause for this difference may be the different sample 
selection criteria, but the trend is not very strong in any of the stud¬ 
ies considered here. 

The same is true for the relation between the depolarisation 
measure ratio and the flux ratio (see Figure lTzl . Neglecting the five 
sources with large flux ratios, for 13 of the remaining 21 sources the 
brighter lobe is less depolarised. Although this is consistent with 
the beaming scenario, the trend is not very significant. 

In summary the trends we find in our sample favour relativistic 
beaming of the hotspots as the underlying mechanism giving rise 


to the observed source asymmetries. However, none of the trends is 
very significant and so our results should be treated with caution. 


5 DEPOLARISATION BEYOND BURN’S LAW 

Depolarisation of radio synchrotron emission by a Faraday screen 
external to the radiation source arises from the inhomogeneity of 
the magnetic field in the screen combined with limited spatial res¬ 
olution of the observations. Let the magnetic field and gas density 
in the Faraday screen be coherent in ‘cells’ of a typical size sq. The 
physical size of the telescope beam at the position of the source 
is t. For t <C so we resolve the Faraday screen and no depolarisa¬ 
tion takes place. In fact, by mapping the rotation measure, we are 
mapping the structure of the Faraday screen in this case. However, 
usually t is comparable to or larger than so and so one beam con¬ 
tains several cells with different position angles of the polarisation 
vector, p. The average polarisation, < |p| >= p, is therefore re¬ 
duced by the partial cancellation of polarisation vectors associated 
with individual cells within the beam. Clearly, the physical size of 
the telescope beam depends on the redshift of a source and so p 
depends in general on z even if the structure of the Faraday screen 
does not change as a function of redshift. 

In previous sections we have used Burn's law (B66) to cor¬ 
rect our measured depolarisation for redshift effects. However, this 
model of the external Faraday screen only applies if the structure of 
the Faraday screen is Gaussian and, more importantly, if so <C 1, i.e. 
for cases where the structure of the Faraday screen is completely 
unresolved. For Cygnus A ICarilli & Barthe! ( 1996 ) report coher¬ 
ent RM fluctuations on scales of jo ~ 10 kpc from observations re¬ 
solving individual cells. At least in this source we have so ~ t and 
Burn’s law does not apply. 


5.1 Alternatives to Burn’s law 

If we cannot assume that so U we need to take into account the 
autocorrelation or structure function of the rotation measure cells in 
the external Faraday screen to derive estimates for the observed RM 
and polarisation properties of our sources. We follow the approach 
bv lTribbl5Nl99ll . hereafter T91), assuming a Gaussian random RM 
field with dispersion a and three different structure functions, D(s), 
where s measures physical distance projected on the plane of the 
sky. Also, because of the limit set by the resolution of our observa¬ 
tions at wavelength X, the RM dispersion we measure, Grmz, is not 
equal to the underlying, ‘true’ dispersion, a. In detail we have: 

• Quadratic structure function with 

D(s) = 2a 2 ( s/sq ) 2 , for s < sq 


D(s) = 

2a 2 , forv > so, 

(2) 

leading to 



A 

-n” 

to 

V 

l-exp(-p-y) 

" 1 + V/p 


a 2 

°RMz 

= cr {1 — [1 — exp (—P)] /|3} . 

(3) 

• Gaussian autocorrelation function with 


D(s) = 2cr 

[l-exp(-s 2 /sjj)] , 

(4) 

leading to 
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o / rad m 2 sq / kpc 


Varying so 

Quadratic 
Gaussian 
Power-law, m — 2 
Power-law, m — 4 

530 

500 

490 

460 

0.15(1 +z) 4J 
0 .11(1+z) 4 - 2 
0.03(1+z) 51 
0.08(1+z) 4 - 4 

Varying a 

Quadratic 

5.0(1+z) 4 - 9 

10 . 


Gaussian 

5.1(1 +z) 4 - 7 

10 . 


Power-law, m — 2 

5.1(1+z) 4 - 5 

10 . 


Power-law, m — 4 

5.1(1+z) 4 - 6 

10 . 


Table 11. Summary of best-fitting model parameters for the dispersion of 
the rotation measure, Orm • 



_ a 

°RMz - 1 + 1 /p' 

• Power-law autocorrelation function with 


D(s) = 2a 2 
leading to 


1 + 2s 2 /msoj 


nj 2 


Redshift 

(5) 

Figure 13. Best fitting models for constant dispersion, o, and varying so- 
Symbols show the data for our sample. Solid line: Quadratic structure 
function. Short-dashed line: Gaussian autocorrelation function. Dotted line: 

( 6 ) 

Power-law autocorrelation function with m — 2. Long-dashed line: Power- 
law autocorrelation function with m — 4. Model parameters are summarised 
in Table 1 11 1 


<|P| 2 > = yPexp(p-y), 

J exp (— m$x/2) exp ^yx~ m l 2 ^ dx 

c RMz = o 2 (mP/2) m/2 r(l-m/2,mP/2) exp (wP/2). (7) 

Here we used the definitions p = Sq/ 2f 2 , y = 4o 2 A 4 and 

r(a,jf) = f exp(—dt (8) 

J X 

for the incomplete T-function. The polarised flux can be obtained 
from p 2 =< |P| 2 > for a Gaussian telescope beam. 

The physical size of the telescope beam at the source loca¬ 
tion, t, is given by the resolution of the telescope and the adopted 
cosmological model. The quadratic and Gaussian structure func¬ 
tions then have two free parameters, the typical cell size, so, and 
the dispersion of the Gaussian random field assumed for the RM 
structure, a. The power-law structure function contains a third pa¬ 
rameter, the exponent of the power-law, m. T91 suggests m = 2 
in his analysis, but other values may apply in different sources. 
Since for the power-law model we have p At 4 /'” (T91), we can 
estimate m for sources with measured polarised fluxes at widely 
spaced wavelengths. Fitting this prediction to the observations of 
a small sample of radio gala xies observed at 1.4, 4.8 and 8.4 GHz 
iAkuior & Garringtorlll995 i), we find 1 < m < 4. In the following 
we use m = 4 and m = 2. 


5.2 Dispersion of rotation measure 

For a given set of model parameters so and a we calculate from 
equations J5} and 0 the expected observed RM dispersion, 
Orm, for each source. We compare this model prediction with the 
observed value of Grmz and calculate the % 2 -deviation of the model 
for each source. We then minimise the overall x 2 -value for the en¬ 
tire sample by varying sq and a. In the following we do not aim for 
a formal best-fitting model. We are only interested in the general 
form of models that are consistent with the data. 

Implicitly we assume here that both model parameters do not 
vary with redshift and are the same for all sources. The expres¬ 
sions for Ormz then lead to Orm z ~ 0 for f < so and Orm z ~ CJ for 



Redshift 


Figure 14. Same as Figure^] but for constant ,vo and varying dispersion, 
a. Model parameters are summarised in Table[U] 


t > sq with a very steep rise in between. Since t increases as z in¬ 
creases, the very rapid transition of Grm from zero to a must occur 
at some finite redshift. As Figure^]shows, our data are not consis¬ 
tent with this behaviour. Incidentally, the usually assumed model 
of B66 would predict the same behaviour of Orm as the models 
of T91 with constant a and so discussed here. Thus our data are 
inconsistent with this model. 

Given our data, clearly either .sq or a or both must vary with 
redshift. Therefore we introduce a third model parameter by setting 
either a = const, and so = s z= o(l + or, alternatively, so = const, 
and a = G. = o(l +zY- 

Figures 1131 and show the model results in comparison 
with our observations. Note here that our models cannot fit all 
data points. However, the plots indicate that all of the models are 
roughly consistent with the data and that either the dispersion, <7, 
or the typical cell size, so, or both must be strong functions of red¬ 
shift (see Table ll 1> . None of the models discussed here fits the data 
significantly better than the others. In fact, the model parameters 
are quite similar between the models. Therefore measurements of 
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Redshlft 


Figure 15. Model predictions for the depolarisation between 1.4 GHz and 
4.8 GHz compared with our observations. All models assume a constant 
dispersion c and a varying typical cell size sq. See TableUJJfor details of the 
model parameters. Solid line: Power-law autocorrelation function with m = 
2. Dashed line: Same as solid line, but for m = 4. Dot-dashed line: Gaussian 
autocorrelation function. Dotted line: Quadratic structure function. 



Redshift 


Figure 16. As Figure^] but for constant so and varying a. 

CJrmz alone do not rule out any of the models. However, it is also 
clear that a model varying both parameters a and so with redshift 
will not provide a significantly improved fit. 

5.3 Depolarisation 

From equations J3}, {5} and 0 we can also derive model pre¬ 
dictions for the polarised flux at our two observing frequencies, 

1.4 GHz and 4.8 GHz. The ratio of the two provides us with a model 
estimate of the depolarisation, DM, which we compare with our ob¬ 
servations. We use the best-fitting model parameters of the previous 
Section, summarised in Table II II 

Figure fT5l demonstrates that models with a variation of the typ¬ 
ical cell size, sg, fail to explain the observations. The expected de¬ 
polarisation significantly exceeds the observed values in all cases. 
Although these models provide adequate fits to the observed dis¬ 
persion of the rotation measure, they do not predict the correct be¬ 
haviour of the percentage polarisation as a function of observing 
frequency. Models with a constant cell size sq, but varying disper- 


X, Y 

rxY 

Student-t 

DM Z , z 

0.68 

6.1 

DM Z ,P 

0.49 

3.8 

DM Z ,D 

-0.23 

-1.6 


Table 12. Spearman Rank coefficients for the correlations between depo¬ 
larisation, DM Z and the fundamental parameters. The depolarisation is now 
adjusted using the best-fitting model with a power-law autocorrelation func¬ 
tion with m = 4. Note that the coefficients are given for depolarisation av¬ 
eraged over the entire source. 

sion a provide a better description of the observations (Figure fT6l . 
The depolarisation predicted by all models at high redshifts is still 
somewhat higher than that observed. However, models for which 
the dispersion of the rotation measure is a strong function of red¬ 
shift are clearly preferred compared to those with a varying typical 
cell size. 

We find that our data are reasonably described by a model as¬ 
suming a power-law autocorrelation function with exponent in = 4 
and a dispersion varying with redshift as a = 5.1 (z + l) 4 ' 6 radm~ 2 . 
We now use this model instead of the model proposed by B66 to 
adjust the measured depolarisation of our sources to the common 
redshift z = 1. The adjusted values DM Z show a strong correlation 
with redshift and a weaker correlation with radio luminosity (see 
Table ll2> . If the observed depolarisation of our sources is caused by 
small-scale structure in a Faraday screen local to the sources, then 
we would expect both cjrmz and DM Z to correlate with redshift. In 
Section [3] we found only Grmz to show this trend, but not DM Z . 
The results of this Section suggest that the depolarisation model 
of B66 used in Section^is inadequate in describing our observa¬ 
tions. For more accurate models we indeed find that both ctrmz and 
DM Z correlate with redshift. This result strengthens our conclusion 
in Section EOl that the predominant effect changing the polarisation 
properties of our sources as a function of redshift is the increasing 
disorder in the structure of the magnetic field in the local Faraday 
screen of the sources. This decrease in the order of the magnetic 
fields is reflected in the strong increase in the dispersion of the ro¬ 
tation measure fluctuations in the Faraday screens, cr, with redshift. 
The change with redshift in the ‘true’ dispersion, a, gives rise to 
the behaviour of the observed dispersion, Ormz, and the observed 
depolarisation, DM : . 


6 SUMMARY AND CONCLUSIONS 

In this paper we investigate in detail radio observations of a sample 
of 26 extragalactic radio sources of type FRII. The original ob¬ 
servations and data reduction are discussed in G04. The observed 
spectral index, a, the structure of the rotation measure, RM, and 
depolarisation, DM, of our sources provide us with information on 
the conditions in the Faraday screens local to the sources. 

We first adopt the cosmological redshift, the low frequency 
radio luminosity and the physical size of our sources as the ‘fun¬ 
damental’ parameters. We find that the radio spectral index shows 
no significant correlations with any of the fundamental parameters. 
The same holds for the rotation measure. The latter result is con¬ 
sistent with our suggestion in G04 that most of the rotation of the 
plane of polarisation occurs in our Galaxy. On smaller scales the 
difference of the rotation measure between the two lobes of each 
source, d RM Z , and the dispersion of the rotation measure, <3 rmz> 
both correlate with redshift and with each other. Thus we conclude 
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that the variations of RM on small scales is caused by a Faraday 
screen local to the source. The redshift correlations also imply that 
the properties of these local Faraday screens change as a function 
of cosmic epoch. There is no independent correlation of c IRM Z or 
Crmz with radio luminosity or source size. 

The small-scale structure of the rotation measure gives rise 
to the measured depolarisation of the sources. We initially use the 
model of B66 for the Faraday screens to convert from the mea¬ 
sured DM to DM-, which we would measure if all our sources 
were located at the same redshift z = 1 ■ The redshift dependence 
of d RM : and particularly Ormz should lead to a similar correla¬ 
tion of DM, with z- Flowever, although such a correlation is present 
in our data, it is not significantly stronger than the independent 
correlation of DM Z with radio luminosity. We argue in Section [5] 
that the model of B66 is inadequate for describing our data and 
weakens the correlation between DM, and redshift. We do not find 
an anti-correl ation of depolarisation measure with size. Some pre- 
vious studies_jStrom|[lV73| jStron^ Jager1ll988t P89: lBest et alJ 
1 1 9991 : llshwara^handra^ taU^98l) find such an a nti-correlation 
while another study does not ( lDennett-Thornd.179%1 ). In most cases 
where there is a contradiction between our result and those pre¬ 
sented in the literature, the differences can probably be attributed to 
different telescope resolutions and/or differ ent o bservin g frequ en- 
cies. However, in the case of P89 and toennett-ThorpeH, 11996 1 we 
have no ready explanation for the disagreement. 

We also study the asymmetries between the two lobes of our 
sources and connections between different kinds of asymmetries. 
The results are somewhat inconclusive, but seem to favour rela¬ 
tivistic beaming of the radio emission of the hotspots as the cause 
for the observed asymmetries. 

We initially employ the commonly used model of B66 for the 
Faraday screens affecting the polarisation properties of our sources. 
However, as T91 points out, this ‘standard’ model only applies for 
observations in which the structure of the magnetic field of the 
Faraday screen is not resolved. We therefore use the more sophisti¬ 
cated models by T91 to derive estimates for ctrmz and DM, which 
we then compare with our observations. We find that only models 
for which the dispersion of the distribution of RM in the Faraday 
screens is a strongly increasing function of redshift can success¬ 
fully explain the observations. All our results therefore point to a 
decrease in the order of the magnetic field in the Faraday screens 
with increasing redshift. If the structure of the magnetic field re¬ 
flects the state of the gas in the source environments, then this 
finding implies that radio-loud AGN at high redshifts are located 
in more turbulent environments compared with their low redshift 
counterparts. 
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